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Abstract 

We study the existence of homoclic solutions for reversible Hamiltonian systems 
taking the family of differential equations + au" — -u + /(u, 6) = as a model. 
Here / is an analytic function and a, b real parameters. These equations are im- 
portant in several physical situations such as solitons and in the existence of "finite 
energy" stationary states of partial differential equations. We reduce the problem 
of computing these orbits to that of finding the intersection of the unstable man- 
ifold with a suitable set and then apply it to concrete situations. No assumptions 
of any kind of discrete symmetry is made and the analysis here developed can be 
successfully employed in situations where standard methods fail. 

Key words: reversible hamiltonian systems, homoclinic orbits, saddle-center 
singularity 



1 Introduction 

Homoclinic orbits have attracted the attention of several authors due to their 
important role as a mechanism leading to chaotic dynamics. This phenomenon 
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was first analyzed by Poincarc, and the study of tlie dynamics in tlie neigli- 
borliood of homoclinic orbits was further developed by Birkhoff, Smale and 
Silnikov (see [1]). 

We say that an orbit is homoclinic to a certain critical set p of a dynamical 
system (it could be an equilibrium point or a periodic orbit) if the orbit is 
bi-asymptotic to this set, that is, \imt±oo 0(^) = P- In this work we concentrate 
on the problem of existence of homoclinic orbits to saddle-center equilibrium 
points, in the context of reversible Hamiltonian systems (see below). Such is- 
sue and its implications are discussed in [2], [3], [4], [5] and [6]. The problem of 
finding homoclinic solutions is related to the existence of solitary waves, spe- 
cially in the presence of surface tension [7], elastic structures [8], [9], [10], [11] 
and spatial patterns in phase transition [12] and [13]. Also, homoclinic solu- 
tions are important to proving the existence of stationary finite energy states 
in partial differential equations [14]. 

In general a Hamiltonian system with a saddle-center equilibrium r does not 
possess homoclinic orbits to r. In order for this orbit to exist it is necessary 
and sufficient that the 1-dimensional stable and unstable manifolds to r, de- 
fined on the same 3-dimcnsional energy surface, intersect. As discussed in 
Section 2, the reversibility of the system will alter completely this situation, 
and more interesting cases arise. We use as a model the equation presented 
in the abstract, but our results are more general. In Section 3 we apply the 
ideas developed here to a specific equation by doing the necessary analytical 
and numerical work considerations. In Section 4 we discuss additional appli- 
cations and report on future simulations to explore in more detail the rich 
consequences of the method here presented 



2 DISCUSSION OF THE METHOD 

Consider a two degrees of freedom family of Hamiltonian systems depending 
on two real parameters a and b, {M,uj, H{a, 6)), where M is a four-dimensional 
C°° manifold, u a symplectic form (closed, non degenerate 2- form over M), 
and H{a,b) : M — > i? is the Hamiltonian function. Let X{M) be the set of 
C°° vector fields over M. Given H there exists a vector field Xh G X{M) 
defined by 

uj{Xh, Y) = dH{Y), for all Y e X{M). 

The flow "i! : RxM ^ Mis defined as = Xh(^'(., x)) and, in symplectic 

coordinates {qi,Q2,Pi,P2), the solution is given by the Hamilton equations 
= f andp, = -||. 
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Our Hypotheses are: 



HI. {M,u, H{a,b)) has a saddle-center equihbrium r = , that is, the hn- 
earized field Xh at r has a pair of real eigenvalues and a pair of pure imaginary 
eigenvalues (non hyperbolic equilibrium point). 

H2. [M,uj, H{a,b)) is reversible with respect to Q, that is, there exists Q : 
M — M where Q is an anticanonical {Q*{uj) = —uj) involution {Q~^ = Q) 
with H o Q = H. We call Q the reversibility of the system. 

Theorem 1 For Hamiltonian systems with a reversibility Q and a saddle- 
center equilibrium r, let x be the set of fixed points of Q. If r & x ^"'^ ^^e 
unstable manifold of r intersects x, then there exists an homoclinic orbit to r. 



PROOF. 

Prom H2 and symplectic properties of oj one has: 

i;toQ^Qo'iJ;_,. (1) 

Let ^ be a solution to the Hamiltonian system such that ^(0) G X- Equation 
(1) implies that ^{t) = ° C(0) = AoQo ^(0) = Q o ^^.^ o ^(0) = Qo ^{-t). 
Since Q(r) = r, if hmt__oo ^(t) = r then hm^^oo ^(t) = limt__oo Q o ^(t) = r. 



Thus the problem of finding an homoclinic orbit is replaced by the search of 
intersection of unstable orbits (in general one-dimensional) with the set x (iu 
general two- dimensional). 

The majority of Hamiltonian systems that possess homoclinic orbits to saddle- 
center equilibrium points also exhibit some kind of discrete symmetry. In this 
case such an orbit is easily found by analysing some Hamiltonian sub-system 
with one degree of freedom. Unfortunately in several physically interesting 
systems this symmetry is unknown or non existent. Our work refers to homo- 
clinic orbits in a class of systems with a kind of reversibility found in important 
physical problems and in differential equations. For such systems there is no 
symmetry that can be used to reduce the number of degrees of freedom and 
the method herein presented is a viable alternative to compute homoclinic 
orbits. 
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3 CONSTRUCTION OF HOMOCLINIC ORBITS IN MODEL 
SYSTEMS 



We use the following fourth order family of differential equations as model 
systems 

+ au" - w + f{u, h) = 0, (2) 
where a e IR, 6 e M, / analytic. We transform 2 into an equivalent system 

u — v 

V— Vy 

(3) 

Pu = -u + f{u,b) 
Pv^ -Pu - av 



with Hamiltonian function H{u, v, Pu, Pv, o, b) = PuV + y + ^ + ^ — F{u, h), 
where F is a primitive of /. 

Equation 2 turns up in several branches of physics, for instance, solitary waves 
in the presence of surface tension [7] . In this case / is approximated by v? and 
parameter a is related to the velocity of the wave. Other important cases in 
which equation 2 is relevant refers to localized patterns in elastic structures 
[9] and spatial patterns in phase transition [13] (in this context 2 is known as 
the Fisher-Kolmogorov stationary extended equation). For more applications 
of solitary waves see [15], [16] and [17]. 

Some authors ([18], [4], [19], [20] and [21]) have looked for the conditions on / 

that guarantee the existence of homoclinic orbits to -u = 0, at least for some 
values of the parameters (a, h) and classes of such functions. However there 
has been no efforts, as yet, for members of such classes, to find the curves 
in parameter space where equation 2 presents homoclinic orbits io u — 
In this context, our interest is in the following problem: given a family of 
functions /, find values (a, h) for which equation 2 has homoclinic solutions 
(j) to the origin, that is hm^^-too 0(0 = (0,0,0,0). For one degree of freedom 
saddle-center hamiltonian system, see [22]. 

One can readily show that system 3 has a reversibility Q : {u,v,Pu,Pv) ' — ^ 
(u, —V, —pu,Pv). whose set of fixed points is x = {{u,v,pu.Pv) \ v = Pu = 0}. 
We define $ as the set given by the intersection of the energy zero level, 
{H = 0}, with a Poincare section defined as {pu = 0}. Here $ is represented 
analytically by the equation ^ + ^ + ^ - F{u, b) = 0. 
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Given b, let p{a) be the intersection of a solution of 3 with $. We expect that, 
given the continuous dependence of solutions with respect to the parameters, 
p{a) intercepts the Poincare section {v = 0} inside $. In this way we determine 
for which parameters the respective solution intersects x, the fixed points of 
the reversibility, a fact that characterizes the homoclinic orbit (figure 1). For 
more properties about homoclinic orbits and their relation with invariant sets, 
see [23], [24] and [25]. 




Fig. 1. Method Illustration 

We illustrate these ideas using a specific equation in the family 2 with a known 
homoclinic orbit to m = 0, that is, u{x) = sech{x). Thus 

u{x) = —sech?{x)senh{x). 

u"{x) = —2sech{x)[—sech'^{x)senh{x)]senh{x) — secK^{x) cosh(a;) = 
= 2sech?{x)senh?{x) — sech(x). 

Resulting in 



u" -u + 2u^ = 0. 



(4) 



Multiplying 4 by u and integrating we obtain the constant of motion 



H(u, u) = — + — = E = cte. 

^ 2 2 2 
« ' 

V{u) 



(5) 



Considering the level curves in Figure 2, we obtain a homoclinic orbit for 
E = 0. 
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Fig. 2. Level Curves for 5 



From 5 we plot the potential V{u) in Figure 3 and get the critical point u = 
for £■ = 0. Classically we obtain motion for V{u) < E = since > 0. This 
fact again guarantees the homoclinic property of the orbit. 




Fig. 3. Potential Function for 5 
Taking the energy level i? = we derive, from 5, 

{uf ~u^ + u^ = 0. (6) 
By 6, the second-order derivative of 4 and a change of coordinates lead to 

a/2 

+ —u" -u+llu^- 12u^ = (7) 
2 

Or, 

+ au" -u + f{u, b) = with /(m, b) = b{llu^ - I2u^) (8) 

which possesses an orbit F given by u{x) = sech{x) homoclinic to u = for 
the homoclinic values (a, 6) = Ij. 
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For model 7, $ = [{u,v,Pu,Pv) if + ^ + T " (t^^ " ^u^)] = o} 
obtain energy zero surface for {a,b) = as shown in Figure 4. The 

structure of the corresponding set which depends only on parameter b, is 
plotted in Figure 5. 




Fig. 4. 3D Energy Zero Surface 
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Evolving system 8 towards the homoclinic point (a, b) 
F in the unstable manifold hits the set x of fixed points of the reversibility, 
"reverting" its behavior and connecting F to the stable manifold, thus charac- 
terizing the homoclinic orbit. For a perturbation of order 10~^ in parameter a 
there will be no intersections and the new orbit does not belong to the stable 
manifold of the equilibrium point, as shown in Figures 6 and 7. 

Using the same procedure developed so far, we build the following equation 
from u{x) = sech^^x): 



- -u" 
4 



u + 3 



65 



-u 



40ii^ 



(9) 
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Fig. 6. 2D View of Homoclinic Orbit (1) and Small Pertubation (2) 
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Fig. 7. 3D View of Homoclinic Orbit (1) and Small Pertubation (2) 
We rewrite 9 as: 



+ au" -u + f{u, b) = where f{u, b) = b { —u^ - ^v? (10) 



From this expression one can show that the orbit F : u[x) = sech?'{x) is 
homochnic to m = for (a, b) = ^, 

We developed an algorithm that runs through all points in a grid of (a, b) 
values, with spacing 10~^, searching for intersections of orbits that belong to 
the instable manifold and the set of fixed points of the reversibility defined 
for the system 10. The result, as shown in figure 8, not only confirms the 
homoclinic value (a, b) = (^—^,3^ as expected, but also find an infinite number 
of additional values with this property. 
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Fig. 8. Systems with Homoclinic Orbits in parameter space 
4 Conclusions 



The computational method developed in this work was based on geometri- 
cal and analytical properties observed in reversible hamiltonian systems. The 
main result contained in the theorem of section 2 does not restrict the hamil- 
tonian dynamics. We have chosen saddle-center equilibriums in view of it rel- 
evance in many applications, but we could extend all results discussed in this 
work to a much wider variety of equilibriums. In [20] we find other examples 
of reversible hamiltonian systems and their applications. 

Applying our method to model 2 with a known homoclinic orbit as seen in 7 
and 10, we already had an indication of its efficiency. Not only the expected 
results were confirmed (figures 6 and 7) but an infinite set on new homoclinic 
values was found and plotted in figure 8, which graphic is similar to the bound 
states distribution in multi-pulse embedded sohtons observed in [26]. This 
phenomenon is known to others authors as "cascade of homoclinic orbits" 
[3], [27] and as "explosion of chaotic sets" [24]. 

The next step is to apply the tool developed herein to others situations of 
interest. We are working with equation + au" — u + f{u,b) = where 
f{u, b) = bv?] this system can be employed as a model of solitary waves in 
presence of superficial tension [7]. We intend to report in a forthcoming work 
the distribution of systems with homoclinic orbits in the space of parameters 
and the corresponding phase transition, as done in [26], [28] and [29], comparing 
our results with other analytical and experimental procedures. We hope to 
provide a "structural" point of view for systems with solitary waves in the set 
of hamiltonian vector fields with two degrees of freedom. 
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